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Abstract 

' In this paper it is shown, that l?-splines are scaling functions for any natural N. 

j . Construction of Haar wavelets, Kotelnikov-Shannon wavelets, nonorthogonal wavelets 

en | based on 5-splines is given. 

+3 ■ 

1. Scaling functions. Let N > 1 - an integer, Z - the set of all integers and L 2 ( 
Hilbert space of square integrable functions. 

Definition 1. Function <f(x) G L 2 (R) is called N-scaling, if it can be represented as 

m! <p(x) = ViV> h n tp(Nx -n), (i; 

in 

<n : 



where coefficients h n , n G Z satisfy to a condition \h n \ < oo. The relationship (1) is 
called the N -scale equation (refinement equation). The set {h n } of coefficients of expansion 
in the equation (1) is called the scaling filter. 

Note 1. If iV-scaling function (p(x) has the compact support of length L, then the sum 
in equation (1) is finite and consist at most L(N — 1) + 1 components. 
The Fourier transform of iV-scale equation is 

b : 

where 



Hn(u) 



Le^"^- ( 3 ) 



N 



Function Hq(u) we shall name frequency function of scaling function ip(x). 
Let's give examples of iV-scaling functions. 

1.1. Haar function. It is defined as characteristic function of an interval [0, 1), 



<po(x) 

Fourier transform of the function is 



i, ie[o,i; 

0, x£[0,l) 
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It is easy to see, that for any natural N 



N-l 1 

Vo (x) = (p (Nx) + (p (Nx-l) + ...(p (Nx-N + l) = Vn V h n (p (Nx-n), h n = —=. 

n=o viV 

Frequency function is 

H ( U ) = 1 (1 + e"- + e"^ + . . . e-**- 1 *) = ^ -^-1) ^^/2 
N y ' N smw/2 

1.2. Kotelnikov-Shannon scaling function. Let's show, that function (p(x) = smnx 
is iV-scaling for any natural N > 1. 

Theorem 1. The function (p(x) = S1 ^ x is N -scale function for any natural N > 1. 
The refinement equation for the function has form: 

i \ / /ATX sin7r(iVx — n) n-^, sin7r(A^x — n) 
¥»(a;) = > <p(n/N) = VN > /i n 5 



tt(Nx — n) ^ ty(Nx — n) 

n£Z v ' neZ v ' 



where 



V N (im\ 
K = sin — . 

nn \N / 



Proof. We shall remind, that function f(x) G L 2 (IR) is called band limited function 
[1], if its Fourier transform f(ou) is equal to zero outside of some frequency band 
The Fourier transform (p{u) of function ip(x) = SVD ^ X has the support on interval [— 7r, 7r] 
and (p(uj) e 1 on this interval [— Then the support belongs to the greater interval 
[—ttN,ttN]. From Kotelnikov theorem [2], band limited function f(x) of a band [— fi, Q] 
can be represented as 

nez 

In particular, for our function (p(x) = if = ttN, then 



/ \ f n \ sin(7riVa; — nn) v-^ ( n \ simr(Nx — n) r-v^, 
<p(x)=y <p \—n ) —\- T r 1 =) - 1 = VN) h n (p(Nx-n). 

nGZ v y nGZ v y nGZ 

1.3. 5-spline scaling functions. Well-known [1], [2], that 5-splines are 2-scaling 
functions. Let's prove, that 5-splines are iV-scaling for any natural N > 1. The 5-spline 
of order zero is Haar function (po(x) (characteristic function of an interval [0, 1)). 5-spline 
<fi(x) of order one with support on interval [0, 2] can be presented [2] as convolution ipo(x) 
with itself, <fi(x) = <po(x) * (po(x). The higher-order _B-splines are defined inductively, 
<£n(x) = <Pn-i(x) * (po(x). Further shall consider, that support of 5-spline (p n (x) is on the 
interval [—(n + l)/2, (n + l)/2] - in case of odd n, and in case of even n - on the interval 
[— n/2,n/2 + 1]. Then the Fourier transform of _B-spline (p n (x) has form [1], [2]: 

where K — - in case of odd n and K — 1 - in case of even n. 
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Theorem 2. B-spline <p n (x) of order n is N-scaling function for any natural N > 1 
and any natural n. The right part of scaling equation (p(x) = v^VEjt h n (p(Nx — k) is the 
finite sum. Filter coefficients hk are find from the formula 

hk= 7=— y] — j-^ r, (7) 

a 

where in case of odd n summation is made on all multi-indexes a = (a , an, . . . a^-i), 
satisfying two conditions: 

\a\ = a>o + cti + • • • + ckjv-i = n + 1 

e*i + 2a 2 + ■ ■ ■ + (N - 1)^-1 — (JV — l)(n + l)/2 = fc, 

£/ie mc?ex varies from k = — (JV — l)(n + l)/2 up to k = (N — \){n + l)/2. 

In case of even n, filter coefficients hk are find from the same formula (7), where sum- 
mation is made on all multi-indexes a = (oto,ai, . . . ajv-i), satisfying two conditions: 

\a\ — ao + ot\ + • • • + ctjv-i — n + 1 

ai + 2a 2 H 1- (AT — l)ajv_i — (JV — l)n/2 = fc, 

£/ie mc?ex A; varies from k = — (JV — l)n/2 up to k = (JV — l)n/2 + 1. 

Proof. First shall consider a case odd n. Then the support of y? n (^) is on the interval 
[— (n + l)/2, (n + l)/2] and the Fourier transform looks like (6), where K — 0. Let's 
substitute (6) in scaling equation (2) and shall find frequency function H (cu), 



(p n {Nu) fsm{Nuj/2)\ n+1 ( uo/2 \ n+l ( 1 sin(^/2) x ™ +1 



(p n (u) V / V. sin (^/ 2 ) / sin(o;/2) 

i „iNu>/2 c -iNu/2\ n+l /i \ n+1 / -, -iNu;\™+ 1 

1 ' < ^ ' 1 i(AT-l)a>/2 » ' 



1 i(JV-l)(n+l)w/2 fi , _-iw _i i_ p -i(JV-l)i^ "+ 1 



-e 



(1 + e"*" + • • • + e-^- 1 )")' 



Since n is odd, the multiplier e M-N-i)(»i+i)/2 [ s an integer degree of e %w . Expanding a degree 

(l + e~ tuJ + h e - l ( Ar - 1 ) w ^ n+1 ) we obtain frequency function H (uj) as the finite sum on 

degrees of a variable z = e~ w , i.e. H (uj) = V N^ 1 h k e~ tku) . Factors of expansion of 
function H (u) give us required filter coefficients hk of a scaling equation. Thus summation 
is made from k = — (JV — l)(n+ l)/2 up to k = (JV — l)(n + l)/2. We shall calculate these 
coefficients. For this purpose we use known analog of a Newton binomial formula 



(z + Zl + ... + Zp ) m = V , 

\a\=m 



'0 • • • 



Then 



(JV-l)(n+l) 

(1 + e -- + . . . + e -i(N-i).y+i = Yl 

k=0 \ a \ =n + l, 

«i+2a 2 + • ■ • +(JV-l)ajv-i = k 



a !o;i! . . . ajy^il 
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H Q (u) = -LjeW-Wn+W* (1 + e -- + • • • + e ~^-^) n+1 = 

(JV-l)(n+l)/2 



fc=-(JV-l)(n+l)/2 |a| = ra + l, U 

ai + 2a 2 + • • ■ + (AT - l)a N -i = k + (N - l)(n + l)/2 

From last expression the statement of the theorem follows. In the second case, when n is 
even, the proof is made similarly. The only difference that support tp n (x) is on interval 
[— n/2, n/2 + 1] and Fourier transform of function ip n (x) has a phase factor e~ lu) l 2 according 
to the formula (6). 

Example 1.3.1. 5-spline ifi(x) of order 1, 



ifi!(x) = 

Fourier transform of the function is 



1 — \x\, x £ [—1, 1], 
0, x£[-l,l] 



^, . /sinc<j/2\ 2 1— cos a; 



\ Uj/2 J LU 
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It is easy to see, that at A" = 3, 



ip^x) = ^(3x + 2) + l(p(3x + 1) + (p(3x) + ^(3x - 1) + ^ V (3x - 2), 

H (u) = 4=— f ie i2u) + 2e lu) + 3 + 2e~ ibJ + e~ i2u ) = - (3 + 4 cos u + 2 cos 2cj) . 
^3^ V ; 9 

At any natural N > 1 the following A^-scaling equation has form 

N - 1 AT I I 

W)= 2^ — ^7 — <PK Nx - n )- 

n=-N+l 

Notes, that the family {(f(x — n), n £ Z} form basis of space V C £ 2 (M), but it is not 
orthogonal. 

Example 1.3.2. Let's consider S-spline <f2(x) of order 2, 



<f 2 (x) 



(x + l) 2 /2, are [-1,0], 

3/4- (x- 1/2) 2 , x £ [0,1], 

(x-2) 2 /2, x£ [1,2], 

0, xi[-l,2\. 



Fourier transform of the function is 

For A/" = 3 filter coefficients of the scaling filter is find from the formula (7). Then, 
1 1 2 7 2 1 1 

^-2 = „ r-i h-i = /lo = - 7=, /ii = - 7=, ^2 = - 7=, ^3 = ~ 7=, ^4 



9V3' 3V3' 3V3' 9V3' 3^' 3^3' 9^3' 
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The scaling relation is 

M%) = IMS* + 2) + §y? 2 (3a; + 1) + §y? 2 (3:r) + 

+ §y? 2 (3:r - 1) + §y? 2 (3:r - 2) + |</? 2 (3x - 3) + ±y? 2 (3:r - 4). 

2. A^multiresolution analysis. Here we will define a basic constructions of multires- 
olution analysis for a case of arbitrary natural coefficient of scaling N > 1, for more details 
about it see [1], [3]. 

Orthogonal iV-multiresolution analysis of L 2 (R) is a sequence of closed nested (tele- 
scoping) subspaces • • • C V_i C V C V\ C V 2 C . . . which has properties: 

1) LijezVj is dense L 2 (R). 

2) n iez ^ = {°}- 

3) f(x) G V} & /(Mr) G V} +1 . 

4) There is a function <^(x) G Vo, called scaling function, that functions (po,n(x) = 
</?(# — n), nfZ form an orthonormal basis of space Vo- 

From properties 3 and 4 implies, that functions 

(p jtn = \/Wip(N j x - n), neZ (8) 

form an orthonormal basis of space Vj for any j. 

From if{x) G V C VI follows that function (p(x) is decomposed on basis of space Vi, i.e. 
it is iV-scaling function, 

(p(x) = VN^2h n (p(Nx -n). 

n 

From the refinement equation follows expansion of base function (8) for any j, k G Z: 

<Pj-l,k(x) = ^ h n-Nk¥j,n( X ) = ^2 k nVj,n+Nk(x). (9) 
n n 

Let's remind [1], that integer translates tp n (x) = ip(x — n) form orthonormal basis of 
subspace V C L 2 (R), if and only if 

\(p(u) + 27m)\ 2 =1 a.e . 

neZ 

Hence [1], [3], if translates <p(x — n) form orthonormal basis of subspace Vo, then frequency 
function H Q (u) satisfy the following property 

N-l 

\Ho(u + 2irm/N)\ 2 = 1 constancy a. e. 

m=0 

2.1. Wavelets. Well-known [1], [3], that in case N = 2 in orthogonal multiresolution 
analysis L 2 (M) each subspace V} of V}+i has orthogonal complement Wj and V, + i = V}©W^-. 
There is a function ip(x), called mother wavelet, that the set of its translates ip(x — n) forms 
an orthonormal basis of space W and functions ipj, n (x) = V2Jip(2 j x — n), n G Z forms an 
orthonormal basis of space Wj for any j G Z. 

In case of arbitrary integer coefficient of scaling N > 1 situation some differs. To scaling 
function (fi(x) corresponds N — 1 wavelets ^(x), . . . ,ip N ~ 1 (x) [1], [3], 

rj,\ x ) = VNj^^MNx - n), k = 1, 2, . . . , N - 1. 

n 
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In frequency space we have: 

= G k (£) (p (^) , fc = l,2,...,JV-l, 

Where Gkiyj) = \J N~ x J2 n g^e~ mu there are frequency functions corresponding to wavelets 
ip k (x). Then the following matrix is unitary [1], [3]: 



/ H (u) H (uj + 2ir/N) ... H (cu + 2rt(N - l)/N) \ 
Gi(J) G^uj + 2-k/N) ... G 1 (uj + 2tt{N- 1)/N) 

\ Gat-iH G N - 1 (oj + 2tt/N) ... G N - 1 (uj + 2tt(N - 1)/N) J 



(10) 



Then we get N — 1 spaces of wavelets: W k , k — 1,2,..., N— 1 and the orthogonal expansion 
for any j 

v j-i = ^®w?e...wf _1 . 

Functions 

rf n (x) = \fm^ k {x){N j x -n), neZ. 

form an orthonormal basis of spaces of wavelets Wf, i = 1, 2, . . . , N - 1. 

2.2. Examples of wavelets. We shall present two examples of iV-scaling functions 
and corresponding wavelets. 

2.2.1. Haar wavelets. For Haar function (p(x) (characteristic function of an interval 
[0, 1)), iV-scaling equation is indicated in 1.1. Nonzero coefficients of the filter h n there are: 
h = hi = ■ ■ ■ = hjsf-i = 1/y/N. 

We shall find wavelet-functions ip l {x), . . . , ip N ~ l (x). Let ip l (x) e Wq. As Wq C Vi, 
then tj} l {x) it is decomposed on basis ipi jn (x) of V\\ ip l — Ylk^nVi^- As ip l (x) _L Vo, then 
for any n we get: (ip\ <fo,n) = 0. Space Vo is nested to V\, hence functions <fo,n(x) also 
are decomposed on basis (pi, n (x) of V\ according to the formula (9): </?o,n = ^2 p hk<fi,Nn+p- 
More in detail 

1 1 1 

V^O.n — —j=fl,Nn H 7=<Pl,Nn+l + 1 7=<Pl,Nn+N-l- 

v A v A v A 



Then the condition of orthogonality ^ l to V becomes 

/ 1 N ^ \ 

\ fc p=0 J 

As the family {ip ljn (x)} form orthonormal basis, then from last equality we shall get 

° l Nn + C% Nn+l + ' ' ' + c jVn+JV-l =0, V tj, i — 1,2, . . . , N — 1. 

From condition of orthogonality of wavelets, (ip l ,ipi) = 5^ follows equalities X^ c n c n = 
We obtain the following set of equations: 

J C Nn + C 7Vn+l + ' ' ' + C Nn+N-l = 0? * 71, i = 1, 2, . . . , AT — 1, , ^ 

I E„« = <%, *,j = l,2,...,iV-l. ^ ' 

The system has a variety of solutions. The elementary solution corresponding to a minimum 
set of nonzero coefficients is ambiguously. A solution will be any set of A" — 1 orthogonal 
vectors in a plane in R N , specified by the equation x + X\ + • • • + £jv-i = 
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To solving of system we shall construct an orthogonal matrix of dimension N, first row 
of which is a vector h = (h , hi, ... , fojv-i) = (1, 1, . . . , 1)/VN. It is possible to take, for 
example, the following iV-dimensional matrix, N, 



( 1 
1 
1 



1 





1 





1 






\ 



V 1 



-N + l J 



and then to make a normalization of these rows. 

For example, if N = 3 we can to select the elementary solution: Cq = 1, c\ — — 1, c\ = 



and Cq 



i, 4 

these vectors: 



1, 4 



-2, remaining coefficients we shall accept equal to zero. Normalize 



Then 



and 



1 



and c n 



1 

7!' 



>/6' 



V2 



i) 2 (x) = -j^P{?>x) + ^^ 3x - 1) - V2<f(3x - 2) 



2.2.2. Kotelnikov-Shannon wavelets. Let's consider Kotelnikov function 

sin7ra; 



TlX 



1, WG[-7T,7rJ, 
0, ^ [— 7T, 7r] 



which, as shown in the theorem 1, is iV-scaling function. Let's construct iV-multiresolution 
analysis. Space V , generated by ip(x — n), n G Z, consist of all square integrable functions 
f(x) for which Fourier transform has support on interval [— ir, 7r]. Space Vj consist of 
functions whose Fourier transforms have the support on interval [—N^tt, 7V%]. Let's find 
wavelets . . . , ^^{x). First we shall define frequency function H (u), and then shall 

find functions Gkiyj) from unitary condition of a matrix (10). Function H (u) can be found 
from the scaling equation, (p(u) = H (uj/N)(p(uj/N), 



H (uj) 



1, ue[-n/N,n/N], 

0, G [-7T, - vr/iV] U [7T/JV, 7T]. 



Outside of interval [— ir, it] function H (uj) extends periodically For k — 1, 2, . . . , N — 1 we 
shall put 

1, \u\e [kn/N, (k + l)n/N], 
0, for other uo G [— 7r, 7r]. 



Wavelet tp k G W 7 ^ is defined from the formula 

= G k {u/N)(p(uj/N) = 



1, |w| G [&7T, (/c + 1)tt], 
0, for other uj. 
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Then for k — 1, 2, . . . , N — 1 we have, 

i poo i p—kir 1 p(k-\-l)7v 

^j k ( x ) = — $ k (tu)e luJX dco=— $ k (tu)e iuJX dcu + — $ k (tu)e tuJX dcu = 

2tt y_oo 2vr y_ (fe+ i )7r 2tt J klT 

sin(A; + l)7r:r — sin knx 

TTX 

This function we shall name Kotelnikov- Shannon wavelets of scale N. 

3. Wavelet-transform. First we shall assume, that we deal with the orthonormal 
wavelets. Then we have scaling function <p(x) and wavelets ^(x), . . . ,ip N ~ 1 (x). Wavelet- 
expansion of signal A = {a n } is [1], [3]: 



ai 



.k = ^2h n a n+Nk , d\ :k = ^2g l n a n+Nk , i = 1, 2, . . . , N - 1. (12) 



where {/i n , n G Z} and {g^, n G Z} there are filters of wavelets <f(x) and ip l (x). Let's 
write last formulas as convolution. For this purpose we shall enter the following coefficients: 
h* n = hT n , g l * = gi n . Then 

Qi.fc = ^ K a Nk-n, d\ k = ^ g l *a N k-n, i = 1,2, . . . , N — 1. (13) 

n n 

Thus, wavelet-expansion is made by the conjugate filters {h* n } and {g 1 *}, with next iV-adic 
decimation (a choice only elements with numbers Nk). Procedure of expansion can be 
repeated, having applied it to set of coefficients cA\ = {a l fc }. Reconstruction of array 
A = {a n } by wavelet-coefficients is 

N-l 

a n = J2 h n-N k ai, k + 9 l n-Nk d \,k- ( 14 ) 

k i=l k 

Last formula can be written also as convolution, having made inverse iV-adic decimation 
of arrays cA\ = {ai,fc} and cD\ = {d\ k , . . . d^ 1 }, 

~ J ai^, if m = Nk, ~ i J d^ k , if m = Nk, 

fll ' m = \ 0, if m not multiply Nk ' 1 - m = \ 0,' if m not multiply Nk 

Then the formula (14) becomes 

N-l 

a n — h n - m ai tm + y^ y^ 9n-m d l,m- (15) 

3.2. Decomposition and reconstruction in not orthogonal case. Formulas of 
wavelet decomposition and reconstruction (12) - (15) usable only for an orthogonal case. 
They are unadaptable, if expansion Vj = Vj-i © © • • • © Wf~^ is not orthogonal and 
functions 



<p jtn (x) = VNiip(Nx - n), M n (x) = VN^ l (Nx - n), n G Z, i = 1, 2, . . . , N - 1 
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do not form the orthonormal systems. Let decomposition of a signal {a k } is made by some 
(nonorthogonal) niters {h k } and {g l k }, % — 1, 2, . . . , TV — 1, by (12), 



CJi 



k — E h n a Nk _ n , d\ k — ^ glflNk-m i — 1,2, . . . , N — 1 



Let's find other filters and {<?£}, i = l,2,...,7V — 1, which lead to perfect reconstruc- 
tion of the input signal by (15). It is convenient to solve the task in terms of formal power 
series. 

Let X(z) = Yln a n zU i s the series corresponding to the signal {a k }. We shall define 
frequency functions H (z) = J2 n h n z n and H^z) = ^2 n g l n z n , i = 1, 2, . . . , N — 1. Then 
their effect on a signal is defined by multiplication: 

X (z) = H Q (z)X(z), Xi(z) = H t (z)X(z), % = 1, 2, . . . , N - 1. 

At wavelet expansion it is still necessary made sample of elements with numbers Nk. In 
terms of formal power series this procedure is reduced to a choice of elements of power 
series with degrees, multiple N. It can be made as follows. Let p = e l2n ^ N , then it is easy 
to see, that the following sum contains only degrees, multiple z N , 

1 (X(z) + X{pz) + X{p 2 z) + ■■■ + Xip^z)) = *N k z Nk = E *N k (z N ) k = A{z N ). 

n n 

Here we used property 1 + p + p 2 + . . . p N ^ = 0, which correctly for any root of degree N, 
p 1, of unit. Thus, in terms of power series wavelet expansion is 

1 N-l 1 JV-l 

A(z n ) = -J2 H o(p Sz ) X (p Sz )> A l {z N ) = -Y,H l {p s z)X{p s z), i = 1,2,... N-l. 

s=0 s=0 

Reconstruction is made by other filters G.i(z) = J2 n 9nZ n , i — 0,1,2, . . . , N — 1, by formula 

N-l 

a n = ^2^Nk a l,k + ^ y^ y 9n-Nk d l.k- 

k i=l k 

In terms of power series it means the following: 

N-l N-l N-l 



G t (z)A t (z N ) = 1 Gi(z) E H l { P s z)X{p s z) = 

N-l /N-l \ 

„ E ( E Gim^z) ) X(p°z) = X(z). 



s=0 \ i=Q 

Therefore it is enough for perfect reconstruction, that it was fulfilled equalities 

N-l N-l 

Gi{z)Hi{z) = N at s = 0, Gi{z)Hi{p s z) — N at s = 1, 2, . . . N - 1. 

i=0 i=0 

It means, that 



/ G (*) ... Gjv-i(s) 
G (pz) ... G N -i(pz) 



V Go^" 1 ^ ... G^-i^" 1 



/ ffoW #o(/^) ••• Hoip^z) \ 

H x {z) Hxipz) ... H^p^z) 

V #jv-i(*) H N ^(pz) ... H N ^{p N - l z) ) 
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N. 



We have obtained the following fact. 

Theorem 3. If a matrix of decomposition filters 



( 



H(z) 



H (z) 
Hi(z) 



H (pz) 
Hi(pz) 



H Q (p N - l z) 
H.ip^z) 



(16) 



\ifAT_i(z) H N ^(pz) ... H^p^z) J 



is nonsingular, then exists the perfect reconstruction of a signal by filters Gi{z), i 
0, 1, 2, . . . , N — 1, which matrix 



G{z) 



( Go(z) 
GAz) 



G (pz) 
G 1 (pz) 



G (p 
Gi(p 



N-l. 
N-l. 



\ 



(17) 



\G N ^(z) G N ^(pz) ... G N ^ 1 (p N ~ 1 z) J 



is inverse transposed to matrix H(z). 

Note. In orthogonal case H(z) is a unitary matrix, and G(z) is a complex conjugate 
matrix to H(z). 

4. 5-spline wavelets. We shall assume, that ip(x) is I?-spline scaling function. As 
is known, it does not define orthogonal multiresolution analysis. It is known also, that its 
frequency function H (z) is polynomial. We shall find wavelets ■0 1 (x) . . . , ip N ~ 1 (x) and its 
polynomial frequency functions, which lead to iV-channel expansion of a input signal and 
perfect reconstruction of the input signal by dual wavelets ^(x) . . . ,ip N ~ 1 (x). 

4.1. General constructions. Let <^(x) is I?-spline scaling function and Hq(z) is its 
frequency function. Let Hi(u), . . . ,Hn-i(z) is frequency (polynomial) functions of filters 
of expansion. According to the theorem 3, for perfect reconstruction it is necessary, that 
matrix (16), where z = e~ lu) p = e t2w / N was nonsingular. We shall include the factor 1/VN 
in frequency functions, as in the formula (2). The matrix (16) has a special type. It is 
possible become free of this special type of matrix H(z) with used Fourier transform on 
cyclical group Z/NZ = {1, p, p 2 , . . . , p N_1 } [3]. We shall define 



1 

N 



It is easy to check up, that the sum on the right depends from w = z N . 
transform is defined by the formula 



The inverse 



iV-l 



H t (z) = J2^AA^)- 



j=0 



Then the last relation can be represented as: 



/ 



H(z) 



A(z 



1 

Z 



1 

pz 



p N ' 1 z 



\ 



N-l 



p N-l z N-l 



^(N-lf z N-l 



A(z N )R(z). 



(19) 



/ 



Matrix A{z ) is already arbitrary nonsingular matrix with polynomial elements. Now 
specific of the matrix H(z) is in the matrix R(z). Setting matrix A(z N ), we can con- 
struct special type matrix H(z) by (19) and also frequency functions Hi(z), . . . ,Hn_i(z) 
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of wavelets, and hence wavelets ^(x) . . . ,^p N ~ 1 (x) themselves. We shall assume, that 
polynomial frequency function H (z) is defined, then the first row of the matrix A(z N ) is 
known, 

z-'Hoiz). (20) 



N 

Z N =W 

Now we construct remaining rows. As each element of first row of the matrix A(w) is a 
polynomial, the first row a(w) = (A 0t0 (w), . . . , Ag^-iiw)) can be represented as 

a(w) = a + ociw + • • • + otg-iw 9-1 , (21) 

where a , a±, . . . , OL g -\ there are g vectors from C^. 

Theorem 4. Let (p(x) be a B-spline and H (z) - its frequency function. If 

a(w) = a + a\W + • • • + a. g -\W 9 ~ x 

is expansion of the first row a(w) = (Aofl(w), . . . , A 0! n-i(w)), obtained by Fourier trans- 
form (20) of frequency function H (z) , then 

a + ai-\ h = • (1, 1, . . . , 1). 

Proof. We shall use an induction with respect to degree p of spline. If p = 0, then 
£>-spline is Haar function (po(x). Then 

H (z) = ±(l + z + z 2 + ... + z N - 1 ). 

It easy to see that A j(w) = 1/N, therefore ct(w) = a = N^ 1 • (1, 1, . . . , 1). Let's assume, 
that the result has been proved for splines (p p (x) of a degree p, and we shall prove it for B- 
splines of a degree p+1. As is known, <p p (x) = ip p -i(x) *ipo(x). Hence (pp~(u) = (tpo(u)) p+1 . 
Therefore from refinement equation 

£o» = Ho (jf) W (£) 
we get that the frequency function H ^ g (z) of £?-spline <f g (x) of degree g becomes: 

H , p (z) = (HMr 1 = j^(l + z + z> + ... + z N -T l ■ 
Let's present H 0jP (z) as 

p N-l p 

H o, P ( z ) = ^2 zkN Yl ak ^ = Yl zkN ( a W + a kA z + ■■■ + a^N^z"- 1 ). 

k=0 i=0 k=0 

Then A 0> j(w) = j^Y1iz n =w z ~^ Hq, p {z) is sample of elements with degrees, multiply z N in 
polynomials H 0jP (z), z~ 1 Hq p {z), , z~ N+1 H 0jP (z). Therefore 

A),o = EaUo zkN a k,o — °o,o + 01,0^ + • • • + a P ,oz pN , 
A o,i = ELo zkN(1 k,i = oo,i + ai,iz N H h a pA z pN , 



A),AT-1 — Efc=0 zkN a k,N-l — ao,N-i + ai,N-iz N + ■ ■ ■ + a Pt N-\z 



pN 
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Hence, 

a o — ( a 0,0? a 0,l; • • • j °0,JV-l), 

<*1 = (^1,0) • • • ) Qi,at-i), ^22) 

a p = (o P ,o, . . . , a P) jv-i). 

From the induction hypothesis, the sum of these vectors is a vector (1, 1, ... , 1)/N. 

Now let's consider £>-spline ip p+1 (x) of degree p + 1. Its frequency function ifo.p+iX- 2 ) 
building from H 0tP (z) by multiplication to a polynomial H (z), 

H , p+1 (z) = ±H M(1 + z + z 2 + --- + z"' 1 ). (23) 

Factor 1/iV in this expression we shall not take into account yet. As well as earlier the 
Fourier transform A j(w) = J2 z n =w z"^ H 0tP+ i(z) is sample of elements with degrees, 
multiply z N in polynomials H 0>p+1 (z), z~ 1 H 0tP+ i(z), . . . , z~ N+1 H 0;P+1 (z). We shall uncover 
brackets in (23) (without factor l/N) and make such sample separately for each part. 

For the factor z° we use Fourier transform to Hq jP {z). Then A j(w) = J2z N =w z ~ j Ho,p( z ) 
is sample of elements with degrees, multiply z N in polynomials H 0p (z), z~ 1 H 0p ( y z), 
z~ n+1 Hq, p (z). Corresponding vector coefficients coincide with found earlier in (22), a ,o = 

«o,i = &i, . . . , ao,p = tt p . From the induction hypothesis, the sum of these vectors is 
equal (1,1,..., 1)/N. ' 

For the factor z l we use Fourier transform to H 0jP (z)z. Then A j(w) = jj Ylz N =w z ~" '^o,p( z ) z — 
~k 12z N =w Z ~ i+A Ho,p( z ) is sample of elements with degrees, multiply z N in polynomials 
H 0jP (z)z, H 0;P (z), z~ 1 Hq >p (z) z~ N+2 H 0jP (z). Let's consider more in detail sample of 
elements with degrees, multiple z N in the element H 0tP (z)z. We have, 

p N-l p 

H 0iP (z)z = ^2z kN ^2 a k,i zl ■ z = zkN ( a k,o z + a ktl z 2 H h a^^^' 1 + a KN _xZ N ). 



k=0 i=0 k=0 



Then 



A ,o = +a , N -iz N + ai , N - lZ 2N +... +a p _ liN _ lZ P N + a p , N ^z^ N , 

Aq,i = ao,o +^1,0^ +a2fiZ 2N + . . . +a Pi oz pN , 



A),JV-1 — ^0,^-2 +01,^-2^ + a '2,N~2 z2N + • • • + a p,N-2 zPN ■ 

Hence, 



a i,o — ( 0, a ,o, a 0jl , a 0i AT_2 ), 

= ( O 0l AT-l, Ol,o, Ol,JV-2 )j 



ai jP — ( a p _i 5 Ar_i, a Pi o, a P ,N-2 ), 

a ljP+1 = ( a P) jv-i, 0, 0, ). 

The sum of these vectors is the vector (1, 1, ... , 1)/N again. 

Let's continue this procedure. For the factor z N_1 is used Fourier transform to H N -i i0jP+ i(z) = 
Ho iP (z)z N ~ 1 . Then A li0 j(w) = ^2 z n =w z^^^Hq^z) is sample of elements with de- 
grees, multiple z N in polynomials H Qp {z)z N ~ l , H 0tP (z)z N ~~ 2 , . . . , H 0jP (z)z H 0jP (z). In this 
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case vector coefficients are 



<2jV-l,0 
OtN-1,1 




0, a 0)0 ), 

OQ.Af-l; a l,0 )) 



OiN-l,p 
a JV-l,p+l 



( Op- 1,1, Qp-1,2 
( Qp,l, Op,2, 



Op-l,AT-l, Qp,0 ), 
Op,JV-l, 0, ). 



The sum of these vectors is the vector (1, 1, ... , l)/iV again. 

Summarizing all obtained vector coefficients of their iV groups, we get, that their sum 
is a vector (1,1,..., 1). Taking into account the factor 1/N omitted earlier, we obtain, that 
the sum of all vector coefficients is equal (1, 1, ... , 1)/N. 

The theorem is proved for 5-spline with support on interval [0,p +1]. In a case when 
the support of 5-spline is centered, as in section 1.3, the theorem is proved similarly. 

Let's return to expansion (21) of first row a(w) = (A 0t0 (w), . . . , A 0t N-i(w)). It can 
be complemented up to a nonsingular matrix with the help simple procedure which is 
generalization of a similar construction in an orthogonal case [3]. As vectors aij in (21) can 
be linearly dependent we shall consider more general case, 



where cko, cti, • . • , a g -i are g linearly independent vectors from po(w) , pi(w) , . . . ,p g -i(w) 
are some polynomials. 

Theorem 5. Let a = (a , a±, . . . , a 9 _i) are g linearly independent vectors from C N 
and Po(w) , pi(w) , . . . ,p g ~i{w) are some polynomials, any of which not vanish on the unit 
circle w = e lt . Then there is polynomial loop A(w) on group GL(N,C) such, that the first 
row of A(w) is Yli=o a iPi{ w )- Thus the degree of loop A(w) is at most the maximum degree 
of polynomials Pi(w) . 

Proof. We shall use induction with respect to maximum degree k of polynomials Pi{w). 
Without loss of generality it is possible to assume, that the sum Yli=o a iPi( w ) ^ s ordered 
on increase of degrees of polynomials Pi(w). If k = then a is some nonzero row vector, 
and we can put A(w) = A, where A is a nonsingular matrix with first row a - 

Let's assume, that k > and that the result has been proved for all smaller k. We can 
suppose, that oi g -\ ^ 0. Let P is an one-dimensional projection to vector cx g ~\ along a 
subspace formed by the first vectors a , a±, . . . , a g _ 2 - Then 



Let's remind, that the projection operator P has the following properties: P 2 = P, P(l — 
P) = 0, (1 - Pf = 1 - P. Now we shall define 

(5{w) = a(w)(l -P + Pg ^(w)- l P) = 

= (a p (w) + aipi(w) H h a g ^ip g -i(w)) (1 - P + p g ^ 1 (wy 1 P) = 

= P g -i{wY l a P + p (w)a (l - P)+ 

+ p 1 (w)a 1 (l - P) + pi(w)p g -i(wy 1 a 1 P H h p g -i(w)p g -i(wy l a g -iP+ 

+ p s . 1 (tB)a s -i(l-P) = 

= + a p (w) + aipi(w) H h a g ^ 2 P g -2(w). 

Therefore /3(w) is a polynomial on w of degrees, smaller k, as deg(p s _2(u>)) < deg(p g -i(w)) = 
k. By induction hypothesis, exists polynomial loop B(w) of degrees deg(p s _2(w)) such, that 



a(w) = a p (w) + aip^w) H h a g ^ip g ^i(w), 



a g - 1 P = a g - 1 , aiP = 0, i = 0, 1, . . . g - 2. 
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first row B(w) is f3(w). Then, supposing A(w) = B(w)(l — P + p g -i(w)P), from the ob- 
tained expression for f3(w) follows, that first row A(w) is a(w). It complete the proof of 
the Theorem 5. 

4.2. P-spline filters of expansion and reconstruction. As is known, frequency 
function for P-spline scaling function is polynomial and does not define orthogonal multi- 
resolution analysis. Constructions mentioned above allow to build nonorthogonal wavelets 
with polynomial frequency functions and to find dual functions which ensure perfect re- 
construction of a signal. If polynomial frequency function H Q (z) is set, it is possible to be 
considered, that the first row of a matrix A(z N ) is known 



Remaining rows can be constructed as in the proof of the theorem 5. As each element of 
the first row of A(w) is a polynomial, the first row a(w) = (Aq j0 (w), . . . , A ^ N _i(w)) can 
be represented as 

a(w) = «o + ot\w + • • • + a g -iw 9 ~ l , 

where cto, «i, . • . , ol q -\ are g vectors from C^. Let Pi is an one-dimensional projection to 
a vector ol 9 -\ along subspace formed by the first vectors a , a±, . . . , a s _ 2 Then 

(3(w) = a(w)(l - P 1 + w'^Pj = 

= (a + a x w + ■■■ + ag-xw 9 ' 1 ) (1 - Pi + w~ 9+1 Pi) = 
= w~ 9+1 aoPi + a (l-Pi)+ 

+ wV(l - Pi) + w 1 w~ 9+1 aiPi + ■■■ + w 9 - l w~ g+l a g _ l P l + 
+ w 9 - 1 a g . 1 (l-P 1 ) = 



= a + a g -i + aiw + h a g - 2 w 



9-1 



Then 



a(w) = (a + a g -i + a±w H h a g ^ 2 w 9 2 )(1 - Pi + w 9 : Pi). 

Continuing this procedure, we obtain the following expansion: 

9-1 

a(w) = (a + a ± + ■ ■ ■ + a fl _i) - P k + w 9 ~ k P k ), 

k=i 



(24) 



where P k is an one-dimensional projection to a vector a g - k along a subspace formed by 
remaining vectors ctj. In the theorem 4 it is shown, that the vector a + «i + • • • + ot g -i is 
(1, 1, . . . , 1)/N. This vector needs to be complemented up to a nonsingular matrix A of 
the order N. We shall specify two most simple methods: 



A 






1 


. 


. 











1 . 


. 











. 


. 


1 



A 



/I 1 

1 -1 

1 1 



1 / 



1 





1 





1 






\ 



V 1 



1 



1 -N + l J 



In the second case the matrix Aq has orthogonal rows and its possible to make an orthogonal 
by normalization of rows. Thus there are many such matrixes Aq. Each such matrix can 
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be obtained from some one by multiplication at the left on an arbitrary nonsingular matrix 
g of the form 

/ 1 ... \ 

di o 0,1,1 ■ ■ ■ Oi n-1 

9 



\ On-1,0 On-1,1 ■ ■ ■ On-1,N-1 ) 

After a choice of constant matrix A , the matrix A(z N ) is given by formula 

1 9-1 

A(w) = -A J](l - P k + w 9 ~ k P k ), 



(25) 



k=i 



where w = z N . Now the matrix of initial filters H(z) become H(z) = A(z N )R(z) ) where 
R(z) is the complementary matrix indicated in the formula (19). The matrix G(z) of 
reconstruction filters is inverse transposed to a matrix initial filters. As H(z) = A(z N )R(z), 
then 



G\z) = H-\z) = R^A- 1 



N 



/I 


z- 1 


1 




V i 


p -(N-l) z -l 



p -(N-lf z -(N-l) 



A(z N ) 



N\-l 



Last formula enables to find filters of reconstruction (dual wavelets). We shall show it on 
examples of 5-splines of order 1 and 2 with parameter of scaling N = 3. 

Example 4.2.1. Let N = 3. We shall consider 5-spline of order I. As we notes earlier, 
frequency function looks like 



H (z) = -z~ 2 {l + 2z + 3z 2 + 2z 3 + z A ). 
9 



Let's calculate 



Z N = W 

Let's put w = e~ tu) and p = e l2w ^ 3 then if zq = e~ %w l 3 summation is made on the following 
values z G {z , pz , p 2 z }. In the further evaluations we shall omit a factor 1/9. Then we 
shall put, 

m ( z ) = z ~ 2 (l + 2z + 3z 2 + 2z 3 + z 4 ) = z~ 2 + 2z~ x + 3 + 2z + z 2 . 
also we shall calculate a Fourier transform on cyclical group of the third order {l,p, p 2 }, 

4>,oH = 3 m °^) = ^ 9 = 3 ' = \Y1 z ~ lm °(z) = ^(3^ 3 + 6) = z~ 3 + 2, 

z' A =w z :i =w 

A , 2 (w) = \J2 z ^ m o( z ) = \( Qz ~ 3 + 3 ) = 2z ~ 3 + 1 • 

z :i =w 

We decompose the first row on degrees with vector coefficient, 

(A ,oH, A ,i(w), Aq >2 (w)) = (3, 2 + z~ 3 , 1 + 2z~ 3 ) = (3, 2, 1) + (0, 1, 2)w~\ 
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a = (3,2,1), ai = (0,1,2), z~ 3 = w' 1 . 

Let P there is an one-dimensional projection to a vector cci along a i.e., «iP = cti 
a P = 0. Now we shall define 

fi(w) = a(w)(l - P + wP) = (a + aiw _1 )(l - P + iuP) = a + «i. 

It is obvious, that (1 — P + w _1 P) is a inverse matrix for (1 — P + wP), therefore 

a(w) = /3(w)(l -P + w~ l P) = (a + «i)(l - P + w^P), where a + a x = (3, 3, 3). 

We complement this first row «o + ol\ up to a nonsingular matrix Ao for example, unit 
rows (0,1,0) and (0,0,1). Then A(w) = Aq(1 — P + w~ 1 P). We shall find expression 
(1 — P + tu _1 P) of an operator. For this purpose we shall add vectors «o = (3, 2, 1) and 
«i = (0, 1, 2) up to basis C 3 . Let a 2 = (0, 0, 1). In this basis the operator (1 — P + w~ l P) 
becomes 

/ 1 
1 - P + w- l P = [0 w' 1 
\ 1 

Let's find expression in standard basis eo, ei, e 2 of spaces C 3 . For this purpose we shall find 
a transition matrix B from one basis, to another a = Be, e = B _1 a, 



B 



Let P ] ak is the matrix of operator in basis (cto, o?i, 0,2)- The operator matrix in standard 
basis is evaluated with help of the formula P e = B~ 1 P a B, P^ s = {B^ 1 ) k s P 3 ak B q - . Then the 
operator (1 — P + w~ l P) in standard basis has the matrix 

/ 1 2/3(1 -vj- 1 ) 4/3(1 -w- 1 ) 
P e (w) = w~ l -2(1 -w- 1 ) 

\ 1 

It allows to calculate polynomial matrix A(w): 

A(w) = A ((l -P + w- x P)) = A P e (w) = 





2 






1 




\ 











Let's calculate the matrix of frequency functions by the formula H(z) = A(z )R(z), 

7710(2) m (p^) m (p 2 z) \ /3 2 + ^ 3 1 + 2^ 3 \ / 1 1 1 
mi(z) mi(pz) mi(p 2 z) = I 2~ 3 — 2(1 — 2 -3 ) -2 P-2 p 2 2 

777.2(2) m 2 (pz) m 2 (p 2 z) J \ 1 / \ z 2 p 2 2 2 p 4 z 2 

We get, m (z) = z~ 2 + 2z _1 + 3 + 2z + z 2 , 7771(2) = z~ 2 + 2z~ l - 2^ 2 , 7772(2) = z 2 . Therefor 
H (z) = \{z~ 2 + 2z- 1 + 3 + 2z + z 2 ), H x {z) = \ (z~ 2 + 2*" 1 - 2z 2 ) , #2(2) = \z 2 . 
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The matrix G(z) of reconstruction filters is an inversion of matrix H(z), G(z) = R(z) l A(z n ) 1 , 

G (z) G 1 {z) G 2 (z) \ ( 1 z- 1 z- 2 \ ( 1/3 -(l+2z 3 )/3 (l-4z 3 )/3 \ 
G Q {pz) Gi(pz) G 2 (pz) =3 1 p- 1 ^ 1 p- 2 z- 2 z 3 -2(l-r 3 ) 

G (p 2 z) Gi(p 2 z) G 2 (p 2 z) / \ 1 p- 2 z- 1 p' A z- 2 J \ 1 y 

Thus it is necessary to take into account, that in wavelet expansion take part the conjugate 
filters of wavelets. It means, that we must to convert a matrix H(z), with a complex 
conjugate elements. As filter coefficient is real, it is enough to make substitution z = e~ lu) 
onz" 1 . We choose elements of the first row and obtain the following filters of reconstruction, 

G (z) = l, Gi(z) = -2z~ 3 + 3z~ 2 - 1, G 2 (z) = -Az' 3 + Qz~ 2 + I - Qz + 3z 2 . 

Thus, reconstruction filters have the following nonzero elements: 

h = V3, 

gl 3 = -2-\/3, 9U = Sv^, 9o = 

^ = -4^3, gl 2 = 6^, 9% = ^, gi = -6>/3, g 2 = 3^3 . 

Example 4.2.2. Let N = 3. We shall consider 5-spline of order 2. Frequency function 
of (p 2 (x) is 

H (z) = ^z- 2 (l + z + z 2 ) 3 = ^ (z- 2 + 3z~ x + 6 + 7z + 6z 2 + 3z 3 + z 4 ) . 
Let's calculate 

A o,M = ^ J2 z ~ jR ^ z )' w e T ' J = °> !> 2 ■ 

Let's put u> = e - * 1 ^ and p = e l27T ^ 3 , then if zq = e~ lu} ^ 3 summation is made on the following 
values z G {z , pz , p 2 z }. In the further evaluations we shall omit a factor 1/27. Then we 
shall put, 

m (z) = z' 2 + 3z~ l + 6 + 7z + 6z 2 + 3z 3 + z 4 
also we shall calculate a Fourier transform on cyclical group of the third order {l,p, p 2 }, 

4>,oH = ^ m °^) = 6 + 32;3 ' 4),iH = I Yl z ~ lm o( z ) =z~ 3 + 7 + z 3 , 

z :i =w z 3 =w 

A Q , 2 (w) = \Y1 z ~ 2m o( z ) = 3z ~ 3 + 6 • 

Z N =W 

We decompose the first row, 

(A),oH, A 0:1 (w), A , 2 (w)) = (6 + 3iu, w" 1 + 7 + m, 3UT 1 + 6) , 
ct-i = (0,1, 3), «o = (6, 7, 6), ai = (3, 1,0), z 3 = w . 

Then 

a(w) = «_iM) _1 + a + aiiw. 
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Let Pi is an one-dimensional projection to the vector ot\ along ct_i and a>o i.e. a\P\ = a>\ 
and a_iPi = 0, a Pi = 0. Now we shall define 

f3(z) — a(z)(l — Pi + w^Pi) = (a-iw^ 1 + a + aiw)(l — Pi +w^ 1 P 1 ) = + a + oti. 

Therefore 

a(z) = (a-iw^ 1 + ao + «i)(l — Pi + wPi). 
Now let P2 is an one-dimensional projection to a vector c*_i along a + «i- Then 

a(z) = (a_i + a + «i)(l - P 2 + w^ 1 P 2 )(l - Pi + wPi). 

We complement the first row a_i + «o + «i = (9, 9, 9) up to a nonsingular matrix Ao, for 
example, by unit rows (0, 1,0) and (0,0, 1). Then A(w) = A (l-P 2 +w- 1 P 2 )(l-P 1 + wPi). 
Let's find expressions of projective operators in standard basis C 3 . First we shall consider 
basis C 3 , consisting of vectors ct_i, a , oti. In this base operators of projection (l—P 1 +wPi) 
become 

/ 1 \ / w' 1 \ 

1 - Pi + wPi = 1 , 1 - P 2 + w^P, = 10. 

\ w J \ 1 / 

Let's find expression of projective operators in standard basis e , ei, e 2 of spaces C 3 . For this 
purpose we shall find a transition matrix B from one basis, to another a = Be, e = B~ 1 a. 

/ 1 3 \ / 2/9 -1/9 5/9 \ 

P = 6 7 6 , P 1 = -2/3 1/3 -2/3 . 
\ 3 1 / \ 5/9 -1/9 2/9 / 

Let Q J ak is the matrix of operator in basis (ao, ai, a 2 ). The matrix of an operator in 
standard basis is evaluated with help of the formula Q e = B~ 1 Q a B, Q q es = (B^ r ) k s Q 3 ak B q - . 
Then operators 1 — P x + wPi 1 — P 2 + w~ 1 P 2 in standard base have matrixes 

-2/3 + 5/3 w -5/9 + 5/9 w \ j I 2/9 w' 1 - 2/9 2/3 w' 1 - 2/3 \ 

2-2w 5/3- 2/3 w , -2/3 w' 1 + 5/3 -2h;- 1 + 2 . 
-2/3 + 2/3 w -2/9 + 2/9 w 1/ \0 5/9^-5/9 5/3«r 1 -2/3/ 

It allows to calculate the polynomial matrix = A (l — Pi + u>Pi)(l — P 2 + -u; _1 P 2 ): 

/ 6 + 3w (1+w 2 + 7w)/w 3(l + 2w)/w \ 

A{w)=\ 2-2w -(2-7w + 2w 2 )/(3w) 2(-l + w)/w . 
\{-2 + 2w)/3 (-l + w)(-5 + 2w)/(9w) (5-2w)/{3w) J 

We compute the matrix of frequency functions by the formula H(z) = A(z N )R(z). The 
first column gives us frequency functions of wavelets ^(x) and ip 2 (x): 

H (z) = ^ (l + 3z + 6z 2 + 7z 3 + 6z 4 + 3z 5 + z (i ) , 
Hi(z) = (-2-6z + 6z 2 + 7z 3 + 6z 4 -6z 5 -2z 6 ), 
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H 2 (z) = — —r (5 + 15z-6z 2 - 7z 3 - 6 z 4 + 6 z 5 + 2 z 6 ) . 
243z 2 v ' 

The matrix G(z) of filters of reconstruction is an inversion of matrix H(z), G(z) = 
R(z)~ l A{z n )^ 1 . Thus it is necessary to take into account, that in wavelet expansion take 
part the conjugate filters of wavelets. It means, that we must to convert the matrix H(z), 
with a complex conjugate elements. As the filter coefficient is real, it is enough to make 
substitution z = e~ lw on z~ x . We choose elements of the first row and obtain the following 
filters of reconstruction: 

G (z) = — (-2 z~ 5 + 6z- 4 -2z~ 3 + 5 z~ 2 -6z- 1 + 2), 

diz) = z 5 (2 z- 8 - 6 z- 7 + 5 z~ 6 -z~ 5 + 3 z~ 4 - z~ 3 - 10 z~ 2 + 12 z' 1 - 4) , 

G 2 (z) = z 5 (6 z~ 8 - 18 z- 7 + 15 z~ 6 - 15 z~ 2 + 18 z' 1 - 6) . 

Note. The obtained expressions depend on a choice of a matrix A . For example, if to 
take a matrix A with orthogonal rows, 




then 

H x (z) = — ^ (8 + 24z-24z 2 - 28 z 3 -24z 4 + 33z 5 + 11 z 6 ) , 

H 2 (z) = —^—(-14-42z + 24z 2 + 28z 3 + 24z 4 -15z 5 -5z 6 ). 
243z 2 v ' 

Reconstruction filters is 

z b 

G (z) = — (8 z- 8 - 24 Z - 7 + 20 z~ 6 - 19 z~ 5 + 57 z~ 4 - 19 z~ 3 + 20 z~ 2 - 24 z~ l + 8) , 

z 5 

d(z) = — (-2 z- 8 + 6 z- 7 - 5 z~ 6 + z~ b -3 z~ 4 + z- 3 + 10 z~ 2 - 12 z' 1 + 4) , 

z 5 

G 2 (z) = — (-10 z- 8 + 30 z- 7 - 25 z -6 - z~ 5 + 3 z -4 - z~ 3 + 20 z -2 - 24 z' 1 + 8) . 
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